function rest = indexfactorial(alpha)
if min(alpha)<=0
    rest=0;
else 
    rest=1;
    for ii=1:max(size(alpha,1),size(alpha,2))
        rest=rest*factorial(alpha(ii));
    end
end